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The Field Generated by an Arbitrary Current 
Distribution Within a Waveguide 

By J. J. Freeman 

Formulas are derived for the electromagnetic field generated by an arbitrary current 
distribution within waveguides of rectangular, circular, and coaxial cross sections. These 
formulas are obtained by generalizing analogous formulas obtained in a previous paper 2 for 
cavities of the same cross section. As a check on the formulas, it is shown that the field 
from an axially directed dipole within a circular guide reduces to the free space value as the 
radius increases indefinitely. 



I. Introduction 

In pursuing certain investigations, it was found 
necessary to derive explicit expressions for the 
electromagnetic field generated by an arbitrary 
current distribution within a waveguide. Since 
the corresponding problem for a cylindrical 
cavity had already been worked out in a previous 
paper (see footnote 2), henceforth designated as 
Paper I, the results of Paper I were generalized to 
apply to a waveguide of corresponding cross sec- 
tion by the simple expedient of allowing the 
length of the cavity to become infinitely large. 
Although the same formulas may be obtained 
more directly by a different method, 3 the present 
paper derives these formulas using the above 
mentioned limiting process since: (a) no explicit 
general solution of the problem has yet been 
published; (b) although the expedient of first 
using a finite cavity for various computations 
involving electromagnetic fields, and then allow- 
ing the dimensions of the cavity to become 
infinitely large has been widely employed, the 
validity of the limiting process is not entirely 
obvious, and the following derivation throws light 
on this process. 

To avoid repetition it will be assumed that 

1 This paper is a revision of Report CRPL-9-3, issued under the same title 
(Feb. 28, 1947). 

2 J. J Freeman, Theory and design of a cavity attenuator, J. Research 
NBS 40, 235 (1948) RP1868. 

3 J. J. Freeman, The noise spectrum of a diode with a retarding field, J. 
Research NBS 42, 75 (1949) RP1953 (see appendix I). 
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Paper I is available to the reader, and the results 
contained therein will be utilized without further 
comment. MKS units are employed. Although 
it was assumed in Paper I that the medium within 
the cavity was lossless, and that the wave number, 
A* — cuy^e, was real, the results are still valid for 
complex values of k. Accordingly, in the follow- 
ing, we assume that the medium has conductivity 
a, so that 

k 2 =O0 2 IJ.e -\-jjJL0HT. 
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Figure 1. Cylindrical cavity. 

II. The Circular Cylindrical Waveguide 

Consider a cylindrical cavity, of length d (fig. 
1) containing current of volume density /(r'), 
where r' designates the position vector of J, and 



193 



where the time variation, e iwt is omitted. In eq 
15 and 21, Paper I, it was shown that 
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In the above, oV represents an element of 
volume containing J, and r is the position vector 
of the point of observation. The m, n,p, 5, and 
X in eq 1 and 2 have the identifying subscripts 
mn implied but omitted. Also, the even, odd 
designations (J) are omitted for convenience, but 
are implied. For the case of a cavity of circular 
cross-section, of radius a, (fig. 1), 

m e n / — sin nd\ T . . /cos nd\ , rt/ , x . 
m — 7 \ cos no) ^^"Un no) ^^^ 

(3) 

where i u i 2 , and i 3 are unit vectors in the directions 
of increasing rfi, and z, and 
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^wd sJv 2 nm —n 2 J n (v nm ) p n = 2 for n 9* 0. 



(4) 
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5=S„ W ==^ where J' n (u nm )=0. (5) 



Also, 



n = — X J„ (Xr) ( /, Mi + /, ) «/* (5r) i 2 , 

v y \sin nd/ \ — cosnd/r v y ** 
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To generalize eq 1 and 2 to the case of an in- 
finite waveguide, with a conducting plug at z=0, 
we let d— > °° , and evaluate the limit of the infinite 
summation over the index Z as rf— > °° . 

Consider the sum over Z in eq 1 , multiplied by 
the factor 1/d. 
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where 
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g* = k 2 -5 2 . 



The factor 1/Vi enters because it is contained in 
the normalizing factor A 2 nm , (eq 4). Let %=lirld. 
Then 

Q 1 /orf ^h sin {2' sin {2 , 1 ^ sings' s in gg 

since AZ=1. Asd-^>o°, A£-^>0, and the sum be- 
comes in the limit the integral, 
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function, 
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f-.9 2 
Let <7i and # 2 be the two square roots of g 2 =k' 2 —8 2 , 
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which lie in the positive- and negative-imaginary 
parts of the complex plane, respectively. Then 

ft=-0i. (12) 

For frequencies above cut-off, when coVe><> 2 , 
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Since, for z>z\ '? — *" 

vanishes over the positive-imaginary pari of the 
plane, as |£|— >°°, the integral can be evaluated by 
taking the residue around the pole g { (fig. 2). Thus, 
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For z<^z' , 2 _ 2 vanishes over the negative- 
imaginary half -plane, as |{|— >a> ? and the integral 
can be evaluated by taking the residue around the 
pole (fig. 2) g 2 . Thus, for 2<2', 



l 



-2—^2- d£= — 2irj- 



e-f 



2<72 



2.9, 



or 



J 



« x,j£(2-2' 



fz:^ #=? ^ liz "" 1 for 2 < 2 '- ( 17 ) 



.9i 




Q 



z^: 



Figure 3. Contour of integration path used by Sommerfeld. 

If o-=0, instead of integrating along the real 
axis of the s-plane it becomes necessary to indent 
the path at g x and g 2 as shown in figure 'A, obtaining 
the contour that Sommerfeld 4 used in the complex 
integral representation of the Green's function for 
the one-dimensional wave equation. 

From eq 11, one recognizes that Si represents 
two plane waves traveling in the + ^-direction 
(for z^>z'), one having its source at z—z f y and the 
other at z— — z' . Also, since m is a vector lying 
in the plane perpendicular to the ^-direction and 
since only transverse elements of current make a 
contribution to the field in eq 1, the first integral 
in eq 11 represents the contribution due to a 
current element at z=z\ and the second integral 
represents tin 4 contribution of its electrical image, 
situated at z= — z\ Both element and image are 
necessary to insure the vanishing of the tangential 
component of E at 2=0, which corresponds to 
the case of a conducting plug at 2=0. On the 
other hand, if the guide is infinitely long in the 
negative z-direction, then the contribution of the 
image at z= — z f must be omitted, and the second 
integral in eq 11 must be neglected. 

Since the case of a perfectly conducting plug at 
2=0 may be easily derived from the case of an 
infinite guide by superposing the electrical image 
with respect to the plane 2=0, the following is 
restricted to the case of a guide extending to 
infinity in both directions. 

Accordingly, from eq 11 and 17, neglecting the 
contribution from the image, 



(18) 



Substituting eq 18 in oq 1, 

__ -*>/« a a g,/;.m(f')e"."-'"m(r) ^ , 

(z^z'). (19) 



Ficttre 2. Contour of integration path in complex plane. 



* A. Sommerfeld, Deutsche Mathematiker— Vereinigung Jahresbericht, 
21, 309 (1912). 
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The various summations over the index I in 
eq 2 may likewise be extended to the case of a guide 
extending to infinity in the positive 2-direction, 
with a conducting plug at z=0. By omitting 
contributions to the field arising from the electrical 
image with respect to the perfectly conducting 
plane at 2=0 of current elements within the guide, 
the boundary conditions are changed from the 
case of a perfectly conducting plug to the case of a 
guide infinite in both directions. Then, if one 
computes the limiting values of the summations 
by the method of residues and omits the contribu- 
tion from the image one obtains: 
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where, for brevity, 



h i =¥-\ i . 
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If eq 8, 18, 20, 21, 22, and 23, are substituted 
in eq 2, one finally obtains 

F - SJ(r') f^^ ft. 



n(r')n(r) ±jn(r')p(r) Tjp(r')n(r) + 



ff 2 P(r')p(r)le^ h ^ z -^\dr f for z^z'. (24) 



Here h { is that square root of P— X 2 that lies in 
the positive-imaginary half-plane. 

As a check on the derivation, it will be shown 
that waveguide eq 19 and 24 reduce to the correct 



value for the free space field, as the boundaries 
recede indefinitely. 

Consider a current filament of length dl, carry- 
ing current / parallel to the axis of the guide 
and situated at the origin of coordinates. Then 
E TE =0, and eq 24 becomes, 
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since only the symmetric modes for which n = Q 
are excited, arid J is normal to n. 

fJ.pdT' = fJ-hJ (\r')dT'=Idl. 

Since n om =— X J (\r)i l =—-j- J o 0^)i u from eq 6, 
and Pom = J (^r)i 3 , from eq 7, eq 25 becomes 
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Equation 26 may be rearranged into the follow- 
ing more convenient form: 

— Tell °° J (\r\eJ h * 

E ™= 2B (grad div+P) S TO^Tl *" (27) 



Consider the sum 



s=S 



J (\r)e ih ^ 



1 hi<L 2 J\(u 



(28) 



There exists an integer M such that for m^>M, 
u om will be large enough so that J (u om ) may be 
replaced by its asymptotic value 5 , 

Jo («-) « y ^ cos (u om - ?5±i »)• (29; 

Since u om is the m th root of J (v> om ) =0, eq 29 yields, 
for m>M, 

u om =Tr(m+3/4:). (30) 

Accordingly, 

J?K«»)=2/TU <M .,m>Af. (31) 

Equation 28 may then be separated into two sums, 

5 J. A . StrattOD, Electromagnetic theory, p. 359, eq 18 (McGraw-Hill Book 
Co., Tnc, New York, N. Y., 1941). 
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&= 2-J h n 2 72/,. N + O 2-J — „2A "AW. (^2) 
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From eq 30, Am= 



> and since X= — -> 
TV a 



Am=- AX. 

7T 



(33) 



Accordingly, eq 32 may be rewritten as follows: 

m J (\r)e^ z 1 - J {\r)\e ih * 

* ~ & hrf J* (u om ) + 2 , £ h ~ ^ 

X=- (Af+3/4) 

(34) 

Now let a— > a> . As a increases without limit 
the first sum in eq 34 approaches zero, AX— >0; the 
second summation becomes, by definition, an in- 
tegral; and the lower limit of integration becomes 
zero. 



Thus, as a—> oo f 
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Substituting from (eq 23a), 



j hl= _VX 2 -F 
into eq 35, one obtains 
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where that square root is to be taken whose real 
part is positive. 



But 
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is the well-known 6 expansion for —p where 

R= ^fr 2 -\-z 2 . Accordingly, eq 27 becomes 

E TM = (grad div +¥) Q^ ^ i 3 ), (37) 

which is the field due to a current element Idl i 3 
in free space. 7 



J. A. Stratton, Electromagnetic theory, p. ~>7<;, eq 17 (McGraw- 
Hill Book Co., Inc., New York, N. Y., 1JM 1 ) . 

7 J. A. Stratton, Electromagnetic theory, p. 431, eq 44 I McGraw- 
Hill Book Co., Inc., New York, N. Y., 1941). 



III. Rectangular Cross Section 

The same general formulas (eq 1 and 2), hold 
in the case of a cavity of length c, and of rectangu- 
lar cross section, whose x and y dimensions are a 




Figure 4. Cavity of rectangular cross section. 

and b, respectively, with the origin at one corner 
(fig. 4). Here, however, the normalizing factors 
and vector functions corresponding to the rectan- 
gular cross section are as follows: 
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Accordingly, when eq 38 through eq 44 are 
substituted in eq 1 and 2, and eq 18, 20 to 23 are 



Arbitrary Current Distribution Within a Waveguide 



197 



substituted for the infinite summations over the 
index I, one gets for the case of a guide infinite in 
both directions, 
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that lies in the positive-imaginary part of the 
complex plane. 

IV. Coaxial Cross Section 

Similar computations for the case of a coaxial 
cross-section guide of inner radius a and outer 
radius b yield the following formulas: 
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that lies in the positive-imaginary part of the 
plane. 
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The unit vectors ij, i 2 , i 3 are defined in figure 1. 
Here 
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where a nm is a solution of 
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and h is that square root of 

h 2 =¥-a\ 



= 

(56) 



which lies in the positive-imaginar} r part of the 
plane. 

V. Conclusion 

The fields within a waveguide excited by an 
arbitrary current distribution have been calculated 
for the cases of circular, rectangular, and coaxial 
cross section by taking the limit of the correspond- 
ing expressions for cylindrical cavities of corres- 
ponding cross section as the length of the cavities 
becomes infinitely great. 

Washington, April 20, 1949. 



198 



Journal of Research 



